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We construct an (a 2 ) -improved overlap-Dirac operator by designing an improved overlap ker- 
nel, based on the Symanzik improvement program. Field rotation terms are also identified to 
improve off-shell amplitudes for both massless and massive fermions. We check the free disper- 
sion relation and propagator, and show that improved results become to close to the continuum 
ones at low momentum region. We test the effect of improvement on the full-QCD gauge config- 
uration and find that the relativistic dispersion relation is satisfied within a few percent error up to 
m a a 0.5. 
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1. Introduction 

Discretization effect is one of the most significant sources of the systematic error in lattice 
QCD calculations. The improvement of lattice action and operators have therefore been extensively 
studied since the early days of lattice field theory. The most well-known and widely used example is 
the clover fermion action [[I]], which removes the 0(a) error in the Wilson's original lattice fermion 
action. According to the Symanzik's improvement program [g], it adds a dimension-five operator 
to the lattice action to cancel the source of error of 0(a) present in the Wilson fermion action. A 
non-perturbative method to tune the parameter in the action has also been established later [Q]. For 
further improvement, one has to add dimension-six and dimension-seven operators consecutively, 
as discussed in [Qj, for instance. These highly improved lattice actions are not so popular in the 
current lattice QCD simulations, since the action contains many terms with parameters to be tuned. 

One of the reasons for the difficulty of designing highly improved lattice fermion operator 
is that the number of operators to be considered is large because of the explicit violation of the 
chiral symmetry in the Wilson fermion action. Indeed, the 0(a) term appears because of the chiral 
symmetry violation, while the chirally symmetric lattice actions do not have this contribution from 
the beginning as one cannot write down the relevant operator of dimension-five while preserving 
chiral symmetry. The same argument applies at 0(a 2m+l ) in general (for m a positive integer). In 
other words, if one starts the improvement program from chirally symmetric lattice actions, the first 
error one encounters is 0(a 2 ), and once it is removed, the next is 0(a 4 ). Therefore, the effect of 
improvement is much more dramatic than in the case of the improvement of the Wilson fermion. In 
fact, the O (a 2 ) -improvement of the staggered fermion has been worked out and used in numerical 
simulations [|5|]. It uses this property of chirally symmetric lattice fermion action. When used for 
heavy quarks, one can greatly accelerate the convergence to the continuum limit. 

In this work we consider the O (a 2 ) -improvement of the overlap fermion The overlap 
fermion preserves exact chiral symmetry through the Ginsparg -Wilson relation [^]. Although the 
numerical cost is high in the practical use of the overlap fermion, dynamical fermion simulations 
have already been performed by the JLQCD and TWQCD collaborations, from which many in- 
teresting physics results have been obtained thanks to its excellent chiral property (for a recent 
summary, see [^]). 

The improvement can be achieved by two steps, i.e. improvement of the action and the field 
rotation. Since the form of the overlap fermion is largely restricted by the Ginsparg-Wilson re- 
lation, improvement of the lattice action is done by modifying the kernel operator to be used to 
construct the overlap operator. To be explicit, we use the fermion action of Eguchi-Kawamoto [§] 
and Hamber-Wu which is called the D34 action in the convention of [Q]. Once we remove the 
Lorentz-violating discretization effects of 0(a 2 ) by this choice of the kernel operator, remaining 
errors can be removed by field rotations. 

2. Formulation of the improved operator 

The overlap operator in the massive case is defined by 

D oy (m q )=(\-—^)D m + m q , (2.1) 
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where the massless operator D ov is given by 



D - = - (,+ ^m)' X=D ~- P a (2 - 2) 



a 



The parameter p controls the large negative mass of the overlap kernel. The conventional choice 
for the kernel operator is that of the Wilson fermion D w , which is 

Av = £(7mV m - \a^) ~ - V + 0{a 2 ). (2.3) 

Near the continuum limit, it reduces to the continuum Dirac operator H) plus the 0(a) error coming 
from the Wilson term. and A M are first- and second-order covariant lattice derivatives, respec- 
tively. 

Near the continuum limit, the overlap operator with the Wilson kernel becomes 

D ov = ID - ^E> 2 + jZ^ D l + £-2 - f {»>*}) + °(« 3 )- (2-4) 

The (9(a) term can be simply removed by a field rotation proportional to D ov , while the 0(a 2 ) 
terms, especially the third term of right-hand side which violates the Lorentz symmetry, cannot 
be removed. The usual overlap operator thus has an 0(a 2 ) discretization error. To remove the 
Lorentz-violating term, we introduce the improved kernel, which is closer to the continuum limit 
D' w ~ E> + 0(a 3 ). Then, the overlap operator takes a simple form up to 0(a A ) errors: 



D' m = H) - ^B 2 + ^ - + 0(a\ (2.5) 



With this operator we can remove the unwanted terms up to and including the 0(a 3 ) term by a field 
rotation proportional to D ov , and the remaining errors start from 0(a 4 ). 

As an improved kernel which has no 0(a) and 0(a 2 ) errors, we use the D34 action. Massless 
D34 action is defined by 

£> D34 = I> M (1 -ba 2 ^) 7m + £ca 3 A2 . (2.6) 

In order to remove the 0(a ) error at tree level, b = 1/6. The parameter c is an arbitrary parameter 
to control the mass of doublers. We take c = 1/6 in the following. For the free case, this action 
has no 0(a) and 0(a 2 ) error, but it is no longer the case once the gauge interaction is turned on. 
In particular, the 0(a) term may arise as radiative corrections, and one has to add another term to 
cancel it. The explicit form of this action on the lattice is 

2 r 

aD D34 = 48 x , y - - £ (1 - Yn)Un, x S x+ ^ y + (1 + y^t/^S^ 



3 > 



(2.7) 
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We now consider the field rotation to remove the remaining discretization effects. Starting 
from the continuum action, / d 4 x \jjr c (x)(E) + m q )y c (x) with fermion fields \\f c an d Yc one may 
define a rotation 

Yc = &cY Wc = W&c (2.8) 



which produces the action / d x \jr(x)r/ oy (mq)\if(x) corresponding to (2.5). Namely, the rotation 
satisfy the relation D' m (m q ) = Cl c (B + m q )Q. c . So far, the rotation matrices Cl c and Q c are written 
in terms of the continuum operator ID. Note that a field rotation does not affect spectral quantities, 
as far as the Jacobian of the transformation is taken into account. The Jacobian may affect the 
renormalization of the gauge coupling at the quantum level but does not matter at the classical 
level. 

There are several choices of the rotations to identify the continuum Dirac operator as the im- 
proved overlap operator up to neglected higher order terms. Since the higher powers of the overlap 
operator, such as D^ v , in the lattice action is computationally expensive in practical simulations, we 
arrange the field rotation so that they vanish in the lattice action. Our choice of the field rotation is 

With this choice, the massive improved operator takes a simple form 



£>ovK) = ( l-^M{m q ,p) )D' W +M(m q ,p) (2.10) 



1 + \ . It means that one can simply use the conventional overlap oper- 
ator in the numerical simulation except that the kernel is improved. Since the rotation operator is 
proportional to E>, the on-shell quantities are unchanged, and off-shell amplitudes are obtained by 
undoing the rotation. To do so, the lattice version of the rotation is given by 

where Q.i and Q.i are the same as £l c and Cl c up to the 0(a 3 ) terms. The off-shell improved 
propagator is then constructed as Q. L D~ l (m q )£lL = (^ + m q) 1 + 0(a 4 ) l , which does not require 
additional inversion of the overlap operator. 



3. Relations at the tree level 



Here, we compare the improved overlap fermion action with the unimproved one at the tree 
level. We consider the dispersion relation 

E(p) = ^p 2 + m q + 0(a n ), (3.1) 

'We note that this construction of the rotation has an apparent problem that the manifest chiral symmetry of the 
form y^Sf (x,y) + Sf(x,y)y^ = is lost. We will discuss on a modification of the lattice action to satisfy this condition 
in future publications. 
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Figure 2: Effective speed of light for the p = (0,0,0) (left) and p = (2tf/L,0,0) (right). The lattice volume 
L = 16 is assumed; 2n/L ~ 0.39. 



which contains the lattice artifact of 0(a n ). The power n is 2 for the Wilson kernel while it should 
be 4 for the improved kernel. Figure [T] shows E(p) for massless (left) and massive (right) cases. 
We can see that the improved operator certainly gives the dispersion relation close to the continuum 
one. To see more quantitatively, in Figure ^ we show the effective speed of light defined by 

^.m^m., ( , 2) 

p 

for p = (0,0,0) (left panel) and p = (2n/L,0,0) at L = 16 (right panel). The results are shown as 
a function of m q a. These plots imply that improved operator indeed very well reproduces the con- 
tinuum dispersion relation with only a few per cent errors up to m q a ~ 0.5, while the unimproved 
operator shows much larger deviation already very close to m q a = 0. 

We also look at the off-shell amplitude (or the quark propagator) at the tree level. We parame- 
terize the quark propagator Sp(p) as Sp(p) = F\(p)p + F2(p)m q after the appropriate rotation £2^. 
We extract F\{p) and F%{p) through 

F i (?) = 7 ^^tr[ipS F (p)] = l + 0(a n ) (3.3) 
4 p^ 

1 i r+mz 

Hp) = j- ^[m q S F (p)} = 1 + {a n ). (3.4) 

4 m L 
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(ap) 2 (ap) 2 

Figure 3: Left panel shows F\ (p) and right panel shows /^(p) versus {ap) 2 . The direction of momentum is 
p= (1,1,1,1) 

In Figure ||, F\ (p) (left panel) and Fj(p) (right panel) are shown. Since the improved operator has 
no 0(a 2 ) term, the slope of the curve corresponding to the improved action vanishes near (ap) 2 = 0. 
These plots are shown for m q a = 0.5. 

4. Non-perturbative test on a dynamical lattice 

We also test the improved overlap fermion action by calculating the meson dispersion rela- 
tion. We use the gauge configurations including 2+ 1 flavors of dynamical quarks generated by the 
JLQCD and TWQCD collaborations [p. The lattice spacing is about a ~ 0.11 fm, and the lattice 
size is 16 3 x 48. Sea quark masses are m uc ja = 0.015 and m s a = 0.080. 

For the valence quark, we use the improved overlap fermion constructed in this work with 
p = 1 .4. We calculate the dispersion relation of the pseudo-scalar meson at several different valence 
quark masses between 0.050 and 0.800 in the lattice unit. 

The effective speed of light is shown in Figure f|. We observe large statistical fluctuations for 
small valence quark masses, as always happens for the correlators with finite momenta. For larger 
quark mass region, we find that the improved operator indeed gives the value closer to unity. Below 
m q a m 0.5, the deviation of the speed of light from 1 is only a few per cent. 

So far, we use the improved kernel as its original form. However, the 0(a) and 0(a 2 ) errors in 
the kernel operator may appear as radiative corrections. We therefore should tune the parameters 
in the action so that these errors vanish, which is left for future works. Also, we are going to extend 
the formulation so that the improved action produces off-shell amplitudes that are consistent with 
the Ginsparg-Wilson relation. 

This work is supported in part by the Grant-in- Aid of the Ministry of Education (No. 21674002). 
Numerical simulations are performed on IBM System Blue Gene Solution at High Energy Accel- 
erator Research Organization (KEK) under a support of its Large Scale Simulation Program (No. 
09-05). 
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Figure 4: Effective speed of light calculated with the improved and unimproved overlap fermion actions. 
These are calculated from two smallest momentum |p| = 0,2n/L. The results are shown for the overlap 
fermion with the Wilson kernel (circles) and with the improved kernel (squared). Open and filled symbols 
represent the data with a point source and with a smeared source, respectively. 
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